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' Abstract 

A new integration technique is presented for systems of linear partial differential equations 
q | (PDEs) for which syzygies can be formulated that obey conservation laws. These syzygies come 

for free as a by-product of the differential Grobner Basis computation. Compared with the more 
obvious way of integrating a single equation and substituting the result in other equations the 
new technique integrates more than one equation at once and therefore introduces temporarily 
fewer new functions of integration that in addition depend on fewer variables. Especially for 
high order PDE systems in many variables the conventional integration technique may lead to 
an explosion of the number of functions of integration which is avoided with the new method. 
A further benefit is that redundant free functions in the solution are either prevented or that 



O ' their number is at least reduced. 

O 

1 A critical look at conventional integration 

X' 

In this paper a new integration method is introduced that is suitable for the computerized solution 
of systems of linear PDEs that admit syzygies. In the text we will call the integration of single 
exact differential equations, i.e. equations which are total derivatives, the 'conventional' integration 
method (discussed, for example, in JI]). To highlight the difference with the new syzygy based 
integration method we have a closer look at the conventional method first. About notation: To 
distinguish symbolic subscripts from partial derivatives we indicate partial derivatives with a comma, 
for example, d xy ti = e hxy . 

To solve, for example, the system 

= f, xx (1) 

= Xf,y+f, Z (2) 

for f(x,y,z) one would, at first, integrate (JTJ) with 2 new functions of integration g(y, z), h(y, z), 
then substitute 

f = xg + h (3) 
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into (J2J), do a separation with respect to different powers of x to obtain the system 

= g, y 

= g,z+h, y 

= h, z 

and solve that to get the solution 

/ = x(az + b) — ay + c, a,b,c — const. 

The main gain of information on which the overall success was based did happen after the substi- 
tution at the stage of separating (J2J) into 3 equations. The integration of (JH) itself did not provide 
new information. The equation = f, xx is more compact than / = xg + h and equally well 
usable in an ongoing elimination process (Grobner Basis computation). (Similarly, in this sense, 
f(x) = asin(x) + 6cos(x) would not provide new information compared to = /" + / as sin and 
cos are only defined as solutions of this ODE.) The main conclusion is: The integration of a single 
equation does not necessarily imply progress in the solution of a system of PDEs, especially if a 
direct separation does not become possible as the result of substituting a computed function. 
This is the case in the example 

= f, yzz Hex) (4) 

= f, xx +f, z . (=:e 2 ) (5) 

discussed in more detail in the next section. Integration of (J2J to / = gi(x, y) + zg 2 (x, y) + gs(x, z) 
and substitution into (J5J) does not yield a separable equation and is therefore not as straight forward 
to utilize as in the first example. 

There is another problem with the conventional method which seems insignificant at first sight 
but becomes severe for high order PDE systems in many independent variables, for example in the 
application in section El 

Substituting / = gi(x,y) + zg 2 (x,y) + g 3 (x,z) into (0) as done in section (|2.2|) and finding 
the general solution for gi,g 2 ,gs is, strictly speaking, a different problem from finding the general 
solution for / of (@]), (J^J! The general solution for g±,g 2 , as determined in section [2~2l will involve 
among other functions the two essential free functions ge(x), gr(x). From the point of view of the 
original system (jlj), © these are redundant functions as they can be absorbed by g 3 . Redundancy 
is an inherent problem of the conventional integration method which has nothing to do of how 
efficient the remaining system after integration and substitution is solved. In section lathis issue is 
discussed in more detail. 

With the new syzygy based integration the situation is very different. Here the decision whether 
to integrate is based on syzygies, i.e. on relations between equations, like 

= {dl + d z ) ei - d y d 2 z e 2 

in the last example and is not based on the form of a single equation. This extra information content 
coming from the syzygies allows the method to perform useful integrations for systems like P|). (0) 
with an instantly useful result. As will be explained further below, syzygy based integration does 
not only integrate one single equation at a time, but in a sense, it performs an integration which is 
compatible with all the equations involved in the syzygy. (More exactly, it integrates all equations 
= P l at once one time where P l are the components of the conserved current of the conservation 
law of the syzygies.) 
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This restrictive 'compatibility constraint' has the effect that the integral involves fewer new func- 
tions of integration which furthermore depend on fewer variables. Consequently fewer new functions 
have to be computed later on which shortens the computation. Also, fewer redundant functions are 
generated which not only avoids the explosion of the number of intermediately generated functions 
but also simplifies the final solution. These effects are especially important for high order PDEs in 
many variables as explained in section El 

The above distinctions between both integration techniques are not purely academic. Section 
19.21 describes how integrations can be combined with eliminations. To apply integrations early in 
the solution process is not new. This strategy has been pursued by the program Crack for nearly 
2 decades. What is interesting and new is how much more beneficial the syzygy based integration 
proves to be compared with conventional integration. In section 19.21 such a comparison has been 
made. One problem has been solved 3 times with a combination of different modules, including 
elimination and conventional and syzygy based integration. The 3 runs differ only in the priority 
of applying these modules and were compared by their running times as well as the number of 
redundant functions in the final solution. 

About the remainder of the paper 

In section 01 the algorithm is described in general and an overview is provided. 

Using the information content of syzygies in the form of conservation laws seems to be the most 
direct and useful way but it is not the only one possible. In section 0] a variation of the algorithm 
is explained which is based on vanishing curls of syzygies. 

Different aspects of the computation of conservation laws for syzygies are the subject of the following 
section. 

The redundancy problem mentioned above is looked at in detail in section |U1 

Even though conservation laws of syzygies might be known, it may not be advantageous to use 
them if the aim is the exact solution of the original PDE-system. In section [7| examples are given. 

A short description of how syzygies are recorded in section 8 is followed by section El introducing the 
'real-life' application which led to the development of syzygy based integration. In three computer 
runs it is shown that this integration method and elimination can be naturally combined for the 
solution of linear PDE systems. 

In the following section the introductory example is continued and both integration methods are 
compared. 

2 An introductory example 

We continue the above example to explain the basic mechanism of syzygy based integration. A 
more complex example is given in section 

2.1 Treated with the new method 

In applying integrability conditions for PDEs systematically, i.e. in computing a differential Grobner 
basis, identities between equations = e a will result that take the form of differential expressions 
with the e a as dependent variables. 

We consider the simple system ©, ©, i.e. 







/ 




3 



= f, xx +f, z . (=e 2 ) 



Assuming, for example, a total ordering > Q of derivatives that implies d x > Q d z and d y > a d z , a 
differential Grobner Basis computation would first eliminate f, xxyzz through cross-differentiation: 

= e2,yzz —£l,xx = f,yzzz ( = '■ 63) (6) 

then a substitution of f, yzz using e\ yields 

= e 3 - e 1 , z 

and a substitution of e 3 using (JBJ) provides the identity 

= e 2 , yzz -(e 1 , xx +e 1 , z ). (7) 

The choice of ordering does not matter here. Any ordering would have resulted in identity 

In this paper we concentrate ourselves to the integration of syzygies, like (J7J), which either have 

the form of a divergence or can be combined linearly to give a divergence = DiP 1 with suitable 

vector components P l (ek) that are differential expressions in the e^. Only in section we outline a 

variation of this principle to deal with a vanishing curl of syzygies. 

The computation of conservation laws of syzygies has several aspects: how to do it in general, 

why the computation of conservation laws for syzygies is a relatively simple task and how to do it 

in less generality but much faster. In the interest of a compact example we postpone this discussion 

to section [SJ 

There are different ways to write as a divergence. We choose any one with as few as possible 
components (here two: P X ,P Z ). This preference is justified towards the end of this section below 
equation (J2"U|) . The question how conservation laws with fewer components are computed is described 
in section El as well. 

We obtain: 

= -ei, xx + (e 2 , yz -e 1 ), z (8) 
= P x lX +P\z (9) 

( f ixyzz );£ (fyxxyz \z ■ (10) 

In the following we will use the vector P % in two representations, first in terms of e^, in our example 
from the syzygy (jSJ): 

P x = -ei, x , P z = e 2 , yz -e l (II) 

and second the representation of P % in terms of the function /, in our example from the identity 
(USD: 

P f ixyzz i P f ixxyz ■ (^^) 

With P satisfying the conservation law condition © we can write P as a 2-dim. curl 

P x = -Q, z , P z = Q, x (13) 
for some potential Q. Using for P l the representation (J12)) we identify 

Q f ixyz ■ 

The existence of differential expressions in unknowns, say f a , for the potential Q is guaranteed 
because all syzygies and all their consequences like = D^P 1 are satisfied identically for any f a . In 
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the appendix B an algorithm DivInt is given that computes potentials Q^{f a ) in general for an 
arbitrary number of independent variables. 

To do the next step in this example, we are reminded that expressions P % {ej) are linear homo- 
geneous in the ej and that they therefore must be zero, i.e. P x = P z = Q, x = Q, z = 0. This means 
that Q is independent of x, z, giving Q = C\{y) and the new equation 

= Q - ci = f, xyz -ci (=: e 4 ) (14) 

with the new function of integration c\ = c\(y). 

Apart from the integral (|T4^) we also get new syzygies. Having on one hand expressions for P % 
in terms of e^e-i due to equations (fTTJl and on the other hand P l in terms of Q,j from equations 
fl~3|) and Q in terms of e 4 from equation (|T1|) we get two new identities 

= P x + Q, z = -e Ux +e A , z (15) 
= P z - Q, x = e 2 , yz -e\ - e 4 , x • (16) 

As equation e± turns up algebraically in at least one of the new identities, this equation = ex is 
redundant and can be dropped. Redundancy of an original equation due to integration need not 
always be the case but it is the case in this example because at least one of P x and P z happens to 
be algebraic in ex (in this case P z ). Identity (fTBj) has already conservation law form. Substituting 
ex from identity (fT7)|) into (fT5|) preserves this form: 

= (-e 4 ,x ), x + (e2, xy -e 4 ), 2 • (17) 

This completes one syzygy based integration step. Because the new system of equations = e 2 = e 4 
obeys the syzygy ()17j) which has a conservation law form with only 2 components P X ,P Z we can 
start another integration step without having to do a differential reduction or cross differentiation 
step. It turns out there are in total 3 more very similar syzygy integration steps to be performed 
which are summarized in appendix A. After these 3 steps the remaining system to solve consists of 
the 2 equations 

= f, xx +f, z (=e 2 ) (18) 

S 2 

= f, y +y c i - y c 2 - xzc i + zc 2 - xc 3 - c 4 (= e 7 ) (19) 
which satisfy the identity 

= -e 2 , y +e 7 , xx +e 7 , z . (20) 

This is a divergence too but now in three differentiation variables. With three non-vanishing P % 
the condition = D{P l has the solution P l = DjQ 1 ^ with more than one non-vanishing Q %: > and the 
condition = P l = DjQ^ has the solution = R l ^ k ,k with free functions R^ k {x n ) = R^ k \x n ) 
where ^ stands for total antisymmetrization. In three dimensions this introduces one new function 
R(x n ) = R xyz through Q xy = R, z , Q yz = R, x , and Q zx = R, y . By performing a syzygy based 
integration again we would solve the remaining equations (|18j h ()19j ) for one function / but also 
introduce one new unknown function R of all variables and therefore not make real progress. This 
is demonstrated in the first example in section [7| These considerations explain why we try to find 
conservation laws of syzygies with as few as possible non-zero P % . 

We return to our example and decide to integrate = e 7 (i.e. (|19p) conventionally because 

• identity (|20p can not be written as a divergence with only 2 terms and 
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• equation (JT§j) can be integrated conventionally with respect to only one integration variable, 
so we will not introduce redundant functions as discussed in the introduction and in section 

El 

To y-integrate equation (fTDJ) we introduce four new functions di(y), . . . , d 4 (y) through q = di, y and 
one new function d$ = d$(x, z) and obtain 

f = — —d\ + —d 2 + xzdi - zd 2 + xd 3 + d A + d 5 (21) 
o 2 

with the only remaining equation (|18|) now taking the shape 

= d$, xx +d$, z . (22) 

A single equation does not have syzygies and the method can not be applied further. What we 
achieved is the integration of equation (HJ) and the change of equation (JSJ) for 3 independent variables 
into equation (}2"2*|) for 2 variables. 



2.2 The same example in a conventional treatment 

For comparison, we solve the system (J3J), (JSJ) again, this time in the conventional direct way. After 
integrating (HJ) to 

/ = 9i(x,y) + zg 2 (x,y) +g 3 (x,z) (23) 
and substitution of / the equation (0) reads 

= 9i(x, y), xx +zg 2 (x, y), xx +g 3 (x, z), xx +g 2 (x, y) + g 3 (x, z), s . (24) 

In equation (|24j) there is no function that does depend on all variables and each variable does 
occur in at least one function. An algorithm for such 'indirectly separable equations' (ISEs) is 
contained in the package Crack (see |Hj and sub-section 19. 2|) . These equations undergo a series of 
differentiations and divisions (producing a list of divisors) 

• to eliminate all functions of some variable, 

• to do a direct separation with respect to this variable, and 

• to use the same list of divisors now in reverse order as integrating factors to back-integrate 
the equations which resulted from direct separation. 

In the case of equation (|2*4*|) a single ^-differentiation eliminates g 3 and allows a direct z separation 
(as gi,g 2 are independent of z) giving = g 2 (x,y), xxy , = g 1 (x,y), xxy +g 2 (x,y), y and through 
back-integration with respect to y further 

= g2(x,y) ixx +g 4 (x) (25) 
= 5-1 0, y) ,xx +92 (x, y) + g 5 (x) (26) 
= g 3 (x,z), xx +g 3 (x,z) iz -zg 4 (x) - g 5 (x) (27) 
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with new functions of integration 54,^5. Renaming <? 4 = ge(x), xxxx , g 5 = g7{x), xx and integrating 
equations (|25|). gives 

02 = -96(x), xx -xg 8 (y) - g g (y) (28) 

Ai 3 

01 = 0e(z) + -^gs{y) + y^(y) + ^io(y) + gu(y) - 07(a) (29) 

= g3{x,z), xx +g 3 (x,z), z -zg 6 (x), xxxx -g 7 (x), xx (30) 

3 2 

/ = 03(2, ^) + ga(x) + -^gs{y) + ~^g%{y) + xg w (y) +gu(y) - gr(x) 

-z(g 6 (x), xx +xg 8 (y) + g 9 (y)). (31) 

The solution (J3TJ) is identical to (}2"T]) and the remaining condition (J3*U|) is identical to (}2"2"|) if we drop 
the redundant functions g 6 , g 7 which can be absorbed by g 3 and substitute g 8 = —di,g$ = d 2 , #10 = 
d 3 , 0n = d^, g 3 = d 5 . A method to recognize redundancy is described in [T3|. It involves the solution 
of an over-determined system of equations which involves even more effort. 

The introduction of redundant functions #6)07 i n the conventional method was unavoidable 
because after reaching system (|25J) - (|27|) with the task to compute gi, . . . ,g$ the information was 
lost that, strictly speaking, not the most general expressions for g 1 , . . . , #5 need to be computed but 
only the most general expression for / = gi(x, y) + zg%{x^ y) + g 3 (x, z). Setting g Q = g 7 = would 
be a restriction for g 2 and g\ in (|28|1. (|2*nj) but is not a restriction for / in (|31jl. 



3 The algorithm in general 

In our notation x l , i = 1, . . . ,p are the independent variables and f a are the unknown functions 
which do not need to depend on all x % . These functions satisfy equations = e a (x n ,fj) where 
j is a multi- index (standing, for example, for n 2 , i.e. d^d^) and where fj stands for a possible 
dependence on f a and any partial derivatives of f a . Total derivatives appear as Di. Summation is 
performed over identical indices. 

The following description is summarized in the overview underneath. The number(s) at the 
start of each item refer to the line number of the corresponding step in the overview. 

(J32|) . (J33|) : For a given system of differential equations ()32j) the investigation of integrability condi- 
tions (e.g. Grobner basis computation) yields identities (J3l?|) . called syzygies. In these syzygies 
the efc take the role of dependent variables. The program Crack has been used to compute 
syzygies for examples presented in this paper but many other computer algebra programs 
are available (for example, RIF 0, diffalg 0,0, |3] diffgrob2 [B]) although only few generate 
syzygies automatically. 

(J34|) : To find conservation laws of syzygies one either can perform a more expensive but general 
search by using the package ConLaw ^2] or other computer algebra software, or one can 
do a more specialized, less general but faster computation as described in section 15.31 In the 
conservation laws as in the syzygies the dependent variables are the e^. 

In order to introduce as few as possible new functions through a syzygy based integration, 
one aims at conservation laws with as few as possible non-zero P % (see discussion towards the 
end of section l2~Tj) . Possible methods to achieve this are described in section IB~2l 

Most often syzygies are very simple expressions and already have a conservation law form. 
Computing conservation laws is not fully algorithmic but it is argued in section lo~T1 that this 
task is relatively simple for under-determined systems of syzygies. 
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(1351) : If a conservation law for the syzygies is known then the following steps can definitely be 
performed. The question is only whether it is beneficial for the purpose of the computation. 
If one has found a conservation law with only 2 components P l then the integration will 
introduce just one new constant and will always be beneficial. If the conservation law has 3 or 
more components P % then at least one new function of all variables will be introduced. In that 
case, if the purpose of the integration is the solution of the PDE system (132)1 then one would 
have to balance how many functions one can solve for due to the new integrated equations (|3T)J) 
against how many new functions are introduced and possible decide not to continue. Examples 
for syzygy based integrations which are useful from the point of solving PDE-systems and 
others that are not are shown in section If usefulness can not be decided at this stage then 
the integration should be performed and decided afterwards. The computational complexity 
of the integration, i.e. of the algorithm DivInt is very low. 

(|36|) : In the computed conserved currents P l (x,ek) we replace the equation names by their 
expressions (}3"2"|) in terms of x, f a . 

(137)) : The resulting P l (x, f a ) in (J36j) are the input to the algorithm DivInt (given in the appendix 
B) to compute a special solution for the potentials = Q^(x,f a ) satisfying P l = DjQ l K 
Here again ^ stands for anti-symmetrization. DivInt works because the kernel of a diver- 
gence DiP 1 is a curl DjQ l i with Q l i = —Qi l and because = DiP % is satisfied identically in 
all f a and their derivatives. 

(|38|) : Because the syzygies = Q m (x,ek) are linear homogeneous expressions in the e^, therefore 
DiP 1 being a linear homogeneous expression in the Q m is also a linear homogeneous expression 
in the e^. Hence the P % are linear homogeneous expressions in the e a . Consequently, we have 
= P l in the space of solutions of the original equations. 1 

(|39|) : On the other hand, the algorithm DivInt computes expressions satisfying P l = DjQ 1 ^ 
identically and therefore = DjQ^ in the space of solutions of the original equations. The 
general solution of this condition for the is shown in (}3T?|) and is the result of the whole 
computation. Its form depends on the number p of no n- vanishing components P l : for p = 2 
a single constant of integration R is introduced for p > 2 one or more functions R^ k (x) are 
introduced. 

(I40j) : The formal integration of = DjQ^ gives new equations whose right hand sides are abbre- 
viated by eij. 

(|41|) : We are instantly able to formulate syzygies which these new equations = satisfy. 

(142)) . (143)) : If any one of them can be solved for one e m (as indicated in (J12I)) then e m = to can be 
substituted in other syzygies and the original equation = e m (x, f a ) can be deleted as it is a 
consequence of the equations e^, ey in lu(x, e^, ey,j ). 

(j44j): 1. As new syzygies have been generated in (jlTj) there is a chance that anyone of them has 
already a conservation law form, like (|15|) . 

2. The substitution of a redundant equation in step ({12} may also lead to a syzygy in conser- 
vation law form, either in the other newly generated syzygy or in any other syzygies. 

1 When computing a differential Grobner Basis the equations in the final basis are also only differential consequences 
of the initial equations and one would not want to delete them. Here the situation is different. — e m has been 
integrated and can be deleted if e m occurs algebraically in other syzygies. 
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3. Finally, there is always the possibility that the new syzygies combined with other syzygies 
take a conservation law form. This would have to be found out by a computation, for example 



using the program ConLaw. 

Given system: = e k (x, f a ) (32) 

Crack — > Syzygies: = fi m (x, e k ) (33) 

ConLaw^ Cons, law form: = D i P i (x, e k ), (34) 

Is CL useful? If not then stop. (35) 

Conserved current: P i = P l {x, e k )\ ek ^ ek{xJa) = P\x, f a ) (36) 

DivInt -> New potentials: P\x, f a ) = DjQ ij with Q ij = Q [ij] (x, f a ) (37) 

Integration of: = P i = DjQ ij (38) 

{R = const in 2 dim 

r--n, x ( 39 ) 
Rv k , k with R^ k = R^ (x) in >2 dim V ; 

New equation names: = j _ ^ fc ^ | =: (40) 

New syzygies: -> = P l (x, e fc ) - e^^- (41) 

Redundancies? e m = a;(a;, e^, e^^- ) — > (42) 

— substitution of e m = a; in any syzygy (43) 

— deleting equation = e m 

return to the determination of conservation laws for syzygies (44) 



The continuation of the introductory example in appendix A is itemized similar to the description 
above. This allows the reader to go through an example and compare it with the overview step by 
step. 

4 An integration based on curls of syzygies 

The described ansatz of extracting information out of syzygies in order to do integrations is not 
the only possible way. In this section we want to provide a different integration method, this time 
based on vanishing curls of syzygies. We will see that it is even more effective than divergence 
based integration but the required structure of the system of syzygies is more special which is 
the reason why it has not been implemented in Crack. Also, the computation of conservation 
laws for syzygies was implemented so far because computer programs, like ConLaw, are available 
to compute conservation laws and because the existence of conservation laws is a relative weak 
condition for syzygies. The method based on curls is shown in the following overview. 
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Given system: 
Syzygies: 



= e k (x, f ) 
= n m (x,e k ) 



Vanishing curl cond.: Q — DjP 1 ^ with P lJ ' — P^(x, e^), 



Curl free tensor: 
New potentials: 
Integration of: 



to new integral (s): 
New equation names: 



= P^(x : e k )\ ek ^ ek{xJa) = P^(xJ a ) 
pv(x, f a ) = D k Q^ k with Q^ k = Q^ k \x, f a ) 
= pv = D k Q^ k 

{R = const in 3 dim 



O ijk (x f a ) — 

[rv* 1 ^ with R^ kl ^R^ kl \x) in > 3 dim 
f Q ijk (xJ a )-R 



Rv 



ki 



■i 



New syzygies: 
Redundancies? 

return to the 



> P^{x^Cg^ &ijkik 

e m = cu(x, e k , eijk,k ) - ► substitution of e m 

— substitution of e m = u in any syzygy 

— deleting equation = e m 

determination of vanishing curls or divergences for syzygies 



The superficial difference between divergence and curl based integration is that P,Q,R have one 
extra index for the curl based method. This method also needs at least 3 independent variables. 
The following two examples involve each 4 independent variables and allow a closer comparison of 
both methods. 

A typical example: 

For 4 unknown functions a, b, c, d depending on x, y,z,t a system of 6 equations 








dix (^it 



-: e 



xyl j 
yz) i 







b, t -d 



■y 



--■ e. 



is given. It has syzygies 



= 


a, z ~ 


C-ix ( 


■ e yt) , 


= 


e xy,y 


"I - &xz,z 


+ e xttt 


= 


^xy,x 


T &yz,z 


+ e ytt t 


= 


&xz,x 


~~ e yz,y 


+ e zt , t 


= 


&xt,x 


~ e yt,y 


&zt,z 








c -b 

^iy u iz 



--■ e xtj 
--■ e zt ) 



which take the form of a vanishing curl: = DjP 1 ^ for P y = leading to potentials Q^ k 

Q xyz = d, Q txy = c, Q xzt = b, Q ytz = a 
and a single new free function of integration R x v zt = g(x,y,z,t). The resulting integrals are 

a = 9,x, b = g, y , c = g, z , d = g, t . 

A related example for a conservation law syzygy: 
In comparison, the typical example using a conservation law syzygy in 4 independent variables 
would involve 6 unknown functions a, b, c, d, f, g and 4 equations, so a less over-determined system: 

= a, y +b, z +c, t (=:ei), = -a, x +d, z +f, t (-■ e 2 ) 
= -b, x -d, y +g, t (=:e 3 ), = -c, x -f, y -g, z (=: e 4 ). 
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The conservation law = e ljX +e 2 , y +e 3 , z +e 4;i gives P l = and potentials 

Q^ = a, Q xz = b, Q xt = c, Q yz = d, Q yt = f, Q zt = g. 

The resulting integrals are 

a = r, z — s,t , b — u,t ~ r , y , c = s jy —u, z , d = r, x —w, t , / = w, z —s, x , d = u, x —w, y 

with new arbitrary functions r, s, u, w. 

If both methods would be applicable, i.e. if the system of syzygies would provide a vanishing 
divergence and a vanishing curl then one would prefer the curl based integration because it makes 
use of more syzygies. 

The last two examples look very artificial but one could exchange the unknown functions 
a,b,c, . . . by any functionally independent expressions, each of them involving at least one different 
function, and the computations and results would be unchanged. 

The remainder of the paper is concerned with divergence based integration. 

5 How to find conservation laws of syzygies 

In order to find a combination of syzygies that is a divergence one could apply computer algebra 
programs ConLaw as described in [T2], [T3] by regarding the syzygies as the equations and the e a 
as unknown functions. In the following subsections we discuss why computing conservation laws 
of syzygies is simpler than computing conservation laws in general, how one can find conservation 
laws with fewer components than independent variables and how conservation laws for syzygies are 
determined in Crack. 

5.1 Under-determination of syzygies 

If one interprets syzygies as PDEs for unknowns e^, then the original equations = ek(x l ,fj) 
are special solutions of these syzygies where the f a play the role of arbitrary functions in these 
solutions. Because at least one of the f a depends on all variables x % (otherwise the original system 
consists only of ISEs to be treated differently, not by checking integrability conditions), the syzygies 
must be an under-determined PDE-system for the unknowns e^. Computing conservation laws for 
under-determined systems of PDEs is an even more over-determined problem. The conservation 
law conditions have to be satisfied identically in a jet space with coordinates x n , e a and all partial 
derivatives of all e a . The more e a occur in the syzygies the more restrictive are their conservation 
law conditions. Another way to see this is that conditions for integrating factors to give conservation 
laws are obtained by applying the variational derivative (Euler-Lagrange operator) to the product 
of integrating factors and syzygies (see [7]). Because there is one Euler operator for each e a we 
get as many conditions as there are different e a . Finally, the more over-determined a system of 
conditions is, the easier it is to solve. Therefore the subtask of computing conservation laws of 
systems of syzygies is usually not a problem. 

5.2 Choosing between different syzygy conservation laws 

The integration of a syzygy = DiP 1 with two derivatives = D X P X + D y P y is always useful but 
not necessarily the integration of a syzygy with more than 2 derivatives because there is at least 



11 



one new function of integration of all variables (see the example in section EJ). Sometimes there is 
a choice allowing to write a syzygy in different forms, for example 

= ei, x +(e 2 , x ) y + e 3lZ 

can also be written as 

= (ei + e 2 , y ), x +e 3 , z . 

To find out whether a conservation law with fewer derivatives exists one has two options. First, one 
can make an ansatz for the conservation law with fewer derivatives and solves the resulting conditions 
(for example, with the programs ConLawI or ConLaw3). Alternatively, one computes the most 
general conservation law involving arbitrary functions. If a conservation law exists which does not 
contain derivatives DjPi, j = m, . . . ,p then = Dj(CPi), j = 1, . . . ,m — 1 is a conservation 
law with an arbitrary function C = C(x m , . . . , x p ). Reversely, finding a conservation law involving 
an arbitrary function C(x m , . . . , x p ) can be exploited to derive a conservation law involving no 
derivatives with respect to x m , . . . , x p as it is described in [TU] . 

5.3 A faster method to find conservation laws 

Methods described above decide whether a conservation law can be built from syzygies, i.e. whether 
there is one in the differential ideal of the syzygies. Computations to decide this general question 
are potentially much more expensive than the other steps of the syzygy based integration which are 
all very quick. In the program CRACK therefore a different, less general but much faster approach 
is taken. Instead of determining whether a linear combination of syzygies exists that makes up 
a conservation law, the program checks each individual syzygy whether it can be written as a 
divergence. 

This is done by using conventional integration to integrate the syzygy with respect to the first 
variable, say x to obtain P x , then integrating the remainder with respect to the second variable, 
say y to obtain P v and so on. A divergence is obtained when no remainder remains after the last 
variable. To find whether the syzygy can be written as a divergence with only two P % the above 
integration is tried at first with all pairs of two independent variables. For example, in the case of 
syzygy © 

&2iyzz ifilixx ~^~&liz ) 

an x-integration gives P x = —e\, x . The remainder e2, yZ z~ e \iz can not be completely y-integrated 
but z- integrated to P z = e 2 , yz 

6 The redundancy problem 

Redundant functions are unavoidably generated as soon as an equation is conventionally integrated 
with respect to at least two different variables, for example, in the integration of = A, x i )X a to 
= A + g{x l ) + h(x 2 ) where g, h depend in addition on all other independent variables occurring in 
the expression A. If A contains n variables x 1 , . . . ,x n then the arbitrariness of g and of h overlap 
to the extend of one function of x 3 , . . . , x n . In other words, if g and h are computed from further 
equations then there will be one redundant function of n — 2 variables in the solution of the original 
problem. 

Let us work out an estimate of how much redundancy is generated when integrating high order 
equations. If the conventional method integrates 

= A } ( x l)mi ) ... ) (a! r >) m n 
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to 

n mi— 1 

where gry are free functions of all variables apart from x % then any two functions gi a ,gib, a ^ b 
have no overlap as their terms g>j a {x l ) a ) gn > (x l ) b involve different powers of x l . Any other pairs of 
functions g a b,9cd, a ^ c overlap. In total there is an overlap within pairs of functions equivalent 
to 

n— 1 n 

m i x m i ( 45 ) 

j=l j=i+l 

functions of n — 2 variables. In the introductory example the integration of = gave rise to 
1x2 = 2 redundant functions of 3 — 2 = 1 variable and in the 'real-life' application in section |U] 
the integration of = c^ X3X3y2yz for c 4 (t, r, Xi, x 2 , x 3 , y 1; y 2 , 2/3) generates an overlap within pairs of 
functions equivalent to 2x1 + 2x1 + 1x1 = 5 functions of 6 variables and for = c^ XlX2XZXiXiyi y 2y2 
even an equivalent of 21 functions of 6 variables. The overlap of two functions is partially also an 
overlap with other third functions and so on and should not be counted twice when trying to account 
exactly for all the redundancy. But this correction concerns the arbitrariness content equivalent to 
functions of less than n — 2 variables so the above formula fi5|) is a good initial approximation of 
redundancy. Keeping in mind that typically a few hundred such integrations may be necessary, the 
severity of the problem becomes obvious. 

Is the redundancy problem an artifact of the chosen examples? 
If one determines higher order symmetries of PDEs then the symmetry conditions may be linear 
PDEs in, say, 30 independent variables (coordinates in jet space). Usually the general solution of 
this overdetermined linear PDE-system involves constants (corresponding to individual symmetries) 
which means that 30 conventional 'successive layers' of integrations would have to be done, each 
'layer' containing integrations that express a function of n variables through functions in n — 1 
variables. In total at least several hundred integrations may become necessary. From this point of 
view the above mentioned application in section |U] to compute C4 is typical. 

Could redundancy be prevented otherwise? 
Partial differential equations may contain symmetries involving arbitrary functions but if not then 
the general solution of the symmetry conditions contains only constants. In that case choosing a 
strictly lexicographical ordering of derivatives in the elimination process the differential Grobner 
basis will involve ordinary differential equations (ODEs). They may not be in the form of total 
derivatives but at least in case they could be integrated, the redundancy problem would not appear 
as each ODE is integrated with respect to only one independent variable. The drawback is that 
Grobner Basis computations are well known to be computationally much more expensive when 
performed with a lexicographical ordering of variables than when performed using a total degree 
ordering of variables. A total degree ordering will provide shorter equations of lower differential order 
but with mixed derivatives, leading to redundancy with conventional integration. The conclusion is 
that even in the special cases where the general solution of the linear PDE system contains essentially 
only constants, the syzygy based integration is superior allowing to use elimination schemes with 
total degree orderings that are more efficient than schemes using strictly lexicographical ordering 
and still being able to reduce the redundancy problem. 

Does syzygy based integration cure the redundancy problem completely? 
In the course of one syzygy based integration all equations = P l are integrated at once one time. 
If = P l (ej) is equivalent to the whole system = e^, or, like in the introductory example (JU),© 
where successive syzygy based integration integrates the system, then redundancy is avoided. If, 
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on the other hand, only a subsystem of equations = e& is involved in = P l {ej) and the result 
of a syzygy based integration has to be substituted in other equations then redundancy may still 
appear as recorded in table 1 in section 19.21 but to a clearly lesser extend. 

Is there another way to determine redundant functions or constants in order to delete them? 
In computations where each free constant in the solution of an overdetermined PDE-system corre- 
sponds to a symmetry or to a conservation law one is interested to determine and drop redundancy 
in order to get an accurate account of their number. For this purpose a method has been developed 
(see ^3]) but this requires the solution of an overdetermined PDE-system on its own and may 
therefore be expensive. 



7 Cases when a syzygy based integration is not useful 

When applying the new integration method to solve a PDE-systrem it not only matters whether all 
steps are algorithmic but also whether its execution is beneficial. Information contained in syzygies 
is useful if it provides a factorization of differential operators. If they do not factorize (for example, if 
they are of first order) then a syzygy based integration can still be useful if more functions are solved 
for than new functions are introduced. If the divergence D^P 1 contains more than two derivatives, 
i.e. the conserved current P l has more than 2 components, then the integral equations ()39|) contain 
at least one new function R 1 ^ of all variables and we may not gain new information from the 
integration if we can not solve for at least 2 functions. This is demonstrated in the following series 
of 3 examples with successively more functions to be solve for and an increasing usefulness of the 
integration. 

Example: 

When computing the Grobner basis of the two equations 

= f, x +f, y (=:e0 (46) 
= f, z (=:e 2 ) (47) 

for a function / = f(x, y) (and in doing that confirming that they are already a Grobner basis) one 
will generate the identity 

= e 2 , x +e 2 , y -e 1 , z . (48) 

From identifying P x = e 2 from (|4*H|) and the general formula P x = D y Q xy + D Z Q XZ together with 
(gH) we identify Q xy = 0, Q xz = f, Q yz = f. With the new function R xyz = c(x, y, z) substituted 
into the formula Q 1 ^ = R 1 -* ,k the new equations are 

= c, z (49) 
= f-c, x (50) 
= f + c, v . (51) 

After a substitution of / from (jHilj) into (JHT|) they are identical to the original set PBj) . (JUJ), only 
now for a function c instead of /. No progress was made. In contrast, for the following two similar 
examples the integration of syzygies is advantageous. 

Example: 
For the equations 

(52) 
(53) 
(54) 






fix ~^~9iy 


(=■■ ei) 





= fiz 


(=: e 2 ) 





= g, z 


(=: e 3 ) 
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the identity 

= e 2 , x +e 3 , y -e 1 , 
results. Integrated in the above manner it gives 









c, x +g 

1 ■ 



(55) 

(56) 
(57) 
(58) 



leaving only equation (p)%j) for c = c(x, y, z) to be solved, an improvement compared to the original 
system f!52|) - ()54|) . In the next example no equations remain to be solved. 

Example: 
For the equations 









the identity 
leads to 



f iz h-ix 

9ix fiy 



= e Ux +e 2)y +e 3 , 



e 2 ) 
es) 









/ + c, x 

9 + C,y 
h + c, z 



(59) 
(60) 
(61) 

(62) 

(63) 
(64) 
(65) 



with an arbitrary function c = c(x, y, z) and no remaining equation. 

In order to incorporate this method of integration into a general program for solving over- 
determined systems the usefulness of integration has to be judged automatically based on the 
number of derivatives in the divergence and the number of functions solved for. But also other 
adjustments to the whole program have to be made. These are discussed in the following short 
section. 



8 Implementation 

Apart from the implementation of the algorithm DivInt as shown in the appendix B, also changes to 
the package Crack were needed in order to automate syzygy based integrations. When checking 
integrability conditions in a Grobner basis computation the program had to keep track of any 
resulting identities (syzygies). This was done in the following way which is conceptually the same 
as the extended Buchberger algorithm (see, for example, the books P and 0). 

To each equation, for example e 3 in (jUJ), we will assign not only a value, like f, yzzz , but also, 
what we will call a 'history-value' or short 'history', i.e. e2, yzz —ei, xx . This history of an equation 
expresses one equation in terms of other equations, i.e. how it was historically computed doing the 
algebraic or differential Grobner basis computation. At the beginning the history of each equation 
e a is e a itself. Whenever a new equation is computed then not only its value but also its history is 
calculated. For example, when in this example f, yzzz is eliminated from equation (jUJ) using equation 
(£Q) then a new equation = is generated where has the value (as all terms cancel) and has 
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the history value e% — e\, z where e% and e± are replaced by their history values. The history of e\ is 
e 1 whereas the history of e 3 is e 2 , yzz —ei, xx giving for e 4 the history e 2 ,y Z z —£i,xx —£i,z as is shown 
in 0. 

In the next section a substantial application is described which is suitable to demonstrate the 
advantages of the new integration method. 

9 The application that led to the development of the syzygy 
based integration 

9.1 The problem 

A problem introduced to the author by Stephen Anco concerns the computation of all conservation 
laws of the radial SU(2) chiral equation in 2 spatial dimensions where the integrating factors are of 
at most 2 nd order. The equation can be written as a first order system for two 3-component vectors 
j(r,t), k(r,t): 

k, 4 = j, r + j x k 

j,t = (rk) >P /r. 

Equation (|67j) is already in conservation law form: 

(rj), t +(-rk), r =0 

and the only other known conservation law (of energy) has zeroth order integrating factors: 

rk ■ [k, t -j, r -j x k] + j ■ [j, t -(rk), r /r] = (j • j + k • k)) , t + (-rj • k) , r = (68) 

The existence conditions for conservation laws below were generated with the program ConLaw2 
described in jT2]. It generates conditions for 6 integrating factors Qi, . . . ,Q$ (like the multipliers 
rk\, rk 2 , rk^, ji, j 2 , J3 on the left hand side of ()68|)). Each of the Qi is an unknown function of 20 
independent variables t, r, j, k, 1 (= j, r ), m (= k, r ), u (= j, rr ), w (= k, rr ). The system consists of 
18 conditions of the form 

Qiiui Qiiwi f 

and 6 conditions of the form 

= Q3,h ^ + Qsih "i 7 " 2 + Qs'h l zr 2 + Q 3 ,i 2 u 2 r 2 + Q 3 ,j 3 l 3 r 2 + Q 3 ,i 3 u 3 r 2 + Q 3 , kl mir 2 + Q 3 , mi wir 2 
+Q3,k 2 m 2 r 2 + Q 3 ,m 2 w 2 r 2 + Q 3 , k3 m 3 r 2 + Q 3 , m3 w 3 r 2 + Q 3 , r r 2 - Q 6 , n kir 2 - Q 6 , jl mir 3 + Q Ml k\r 
-Q&,h m i 7 ' 2 ~ Qe,h wir 3 - 2Q 6 , Ul ki + 2Q 6 , Ul m x r - Qe, Ul wir 2 - Q 6 , J2 k 2 r 2 - Q e , h m 2 r 3 + Q 6 , h k 2 r 
-Qe,i 2 m 2 r 2 - Q 6 , h w 2 r 3 - 2Q 6lU2 k 2 + 2Q 6 , U2 m 2 r - Q 6 ,„ 2 w 2 r 2 - Q 6 , J3 k 3 r 2 - Q 6 , j3 m 3 r 3 + Q 6 ,i 3 k 3 r 
-Qe,l 3 m 3^ 2 ~ Q6,h w 3 r 3 - 2Q 6 ,„ 3 k 3 + 2Q 6 , U3 m 3 r - Q e ,u 3 w 3 r 2 ~ Q 6 , kl hr 3 - Q 6 , kl j 2 k 3 r 3 + Q 6 , kl j 3 k 2 r 3 
-Qe, mi uir 3 - Qs,mi j2m 3 r 3 ~ Q 6imi hk 3 r 3 + Qe, mi j 3 m 2 r 3 + Q 6 , mi hk 2 r 3 - Qe, Wl j 2 w 3 r 3 - 2Q 6 , Wl l 2 m 3 r 3 
-Q&,wi u 2 k 3 r 3 + Q 6 , Wl j 3 w 2 r 3 + 2Q 6 , Wl l 3 m 2 r 3 + Q 6 , Wl u 3 k 2 r 3 + Q 6lk2 jik 3 r 3 - Q 6lk2 l 2 r 3 - Q 6 , k2 j 3 k x r 3 
+Q6,m 2 jim 3 r 3 + Q e ,m 2 hk 3 r 3 - Q 6 , m2 u 2 r 3 - Q 6 , m2 j 3 mir 3 - Q 6 ,m 2 hhr 3 + Q 6 , W2 jiw 3 r 3 + 2Q 6 , W2 hm 3 r 3 
+Q erW2 uik 3 r 3 - Q 6 , W2 j 3 wir 3 - 2Q e , W2 l 3 mir 3 - Q 6lW2 u 3 kir 3 - Q 6 , fe3 jik 2 r 3 + Qe,k 3 j 2 hr 3 - Q 6 , fe3 l 3 r 3 
-Q6,m 3 jim 2 r 3 - Q e ,m 3 hk 2 r 3 + Q 6 ,m 3 j2 , m 1 r 3 + Q e ,m 3 hfar 3 - Q e ,m 3 u 3 r 3 - Qe, W3 j\w 2 r 3 - 2Q 6 , W3 hm 2 r 3 
-Qe, W3 uik 2 r +Qe, W3 j 2 wir +2Q 6 , W3 l 2 mir +Q§, W3 u 2 k 1 r - Qe.tr -k x Q 2 r +k 2 Qir 
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After introducing new unknown functions x±, %ji through Ui = Xi + jji, Wi = Xi — %ji the 18 short 
equations took the form of a total derivative and each one could be integrated on its own but when 
the computed functions were substituted only indirectly separable equations (ISEs) like were 
obtained. 2 

Despite of the initial success in performing these integrations all attempts to complete the 
solution of the over-determined system failed with the 1999 version of Crack. That this was not 
simply a matter of lacking computing power became obvious after extracting a small sub-system 
of equations for only one of the unknown functions 3 c±{t, r, x%, X2, £3, yi, 2/2, Hz) where some of the 
equations are easy to integrate: 

— c 4 5x3x31/24/3 — C 4;xix'24/i4/34/3 — C 4;xiX24/i4/i yz ~ c Aix\X2Xzy\y\ — c 4iX2Xzxzxzyiyiyz 

— c 4}Xix 2 xzxzx 3 yiy2y2 = c ^ix\X2Vixzy\y\yiyi = ^ix^xzxzyzyzyz ~ c ^)X2Xzxz%zyiyzyz 
C-\ix\X2.xzxzy\yiyi '^' ( -'^ix\x2X2xzy\y2yz ^4ixiX2xzxzyiy2 ^^iix\X2X2xzyiyz ^45x1x2x2x3x31/14/3 %3 

— C4,xix 2 X3X3X34/i2/3 X\ — C4,xix 2 x 3 X3i/3j/3 +£45x2x3x3x34/14/3 

= c 45X3X 3 X3j/ij/3 x l + C 4;X3X3j/ij/34/3 Hi ~ C ^ixzx z yzyz (69) 

= C 4 5X1X2X3X34/11/24/24/2 2/3 + 204,2^ X2X2X24/1J/24/3 _ 2 C 45XiX2X2X 3 J/iJ/24/2 ~!~ C 4!XiX2X2X3X34/i4/2J/2 X 3 

C 4 5X1X2X3X34/14/22/2 V'i 2C4,XiX2X2X 2 X3j/i4/3 ^3 "I" '^' C A:lX\X2X2X2y\yz ^ C ilX l X2X2Xzyiy2 

C4 ?:ria ;2X3X3X34/ll/2 "^1*^3 3C4, Xia ;2X3X34/lJ/2 "^1 ^^AlX\X2X2Xzy\yz "^1 
_ 2 1 2 

C4,XiX2X2X 3 X3i/3j/3 3^3 T C k)X2X2 x Z x 3X3VlD3 X 3 

Even the solution or at least simplification of this sub-system was not possible. The problem was 
not to find equations with the form of a total derivative and to integrate them. The problem was 
the growing number of new functions of integration (which did still depend on 7 variables) and the 
appearance of too many only indirectly separable equations (ISEs). 

Since 1999 the module for handling ISEs has been improved considerably. The current version 
of Crack (Dec. 2001) can simplify the above system quickly using the conventional integration 
of total derivatives. Nevertheless, by adding the ability of performing syzygy based integrations 
the computation speeds up further and the solution involves fewer redundant arbitrary functions. 
Tests described below show that syzygy based integrations are well suited to be performed along 
the computation of a differential Grobner basis without the negative side effect of introducing too 
many redundant functions. By that Grobner basis computations can be cut short and the risk of a 
memory explosion be lowered. 

9.2 A comparison of three computer runs 

Before describing the details of 3 different computer runs, a few comments about the setup have to 
be made. The package Crack for solving and simplifying over-determined PDE-systems contains 
about 30 modules for different actions to be taken either with individual equations or with groups 
of equations of the system. Modules used to solve systems like ()69|) are 

1. Direct separation of an equation with respect to some variable that occurs only explicitly in 
the equation. 

2 Although each of the ISEs is over-determined on its own, this over-determination can not be utilized easily 
because there is no independent variable which occurs only explicitly that would lead to direct separations. 
3 New constants and functions of integration are all called c; in Crack with successively increasing subscript. 
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2. 



Substitution of a function / either by zero or by at most 2 terms and only if other functions 
occurring in these 2 terms depend on fewer variables than /. 



3. Integration of an equation if it consists of a single derivative with respect to only one variable. 

4. Elimination of a function / from any equation if / occurs only algebraically and linearly and if 
/ depends on all variables occurring in this equation. Substitution of / in all other equations. 

5. Deleting of any redundant equations as described on the bottom of the overview in section El 

6. Integration based on a syzygy in conservation law form. 

7. Conventional integration of a PDE but only if sufficiently many integrations are possible such 
that the integrated equation can be used for a substitution. 

8. Indirect separation of an equation (ISE). (This is a complex step which can invoke other direct 
separations and indirect separations of resulting equations.) 

9. Reduction of the leading derivative of one equation with the help of another equation or 
formulation of an integrability condition between two equations. (This is a typical step in a 
Grobner basis computation.) 

10. Any integration of any equation even if not complete. 

These modules are called in a specific sequence which can be chosen by specifying a list of numbers, 
each number representing one module. For example, if in table 1, column 2 the priority list of run 
1 is chosen to be 1 2 3 4 8 9 7 10 then the modules as numbered above are tried in this order until 
one module is successful and then they are again tried beginning with 1 and so on. This is only a 
simplified description of the operation of Crack but it is sufficient for the purpose of this section. 



run 


priority list 


#of 


time 


# of terms 


$ of redundant functions 




of actions 


steps 


in sec 


in equ. (I7T1) 


of 6 var. 


of 5 var. 


of 4 var. 


of 3 var. 


1 


1 2 3 4 8 9 7 10 


1077 


124 


6 


7 


16 


2 





2 


1 2 3 4 7 8 9 10 


1175 


122 


12 


4 


45 


23 


5 


3 


123458697 10 


362 


23 


8 


2 


19 


2 






Table 1. A comparison of three different runs on the system (|69|) . 



In table 1 three computer runs are compared. Column 3 gives the number of successful calls 
of the modules in the priority lists. Times shown in column 4 have been measured in a session 
of the computer algebra system REDUCE version 3.7 with 120 MByte memory (although only a 
few MByte are needed for this computation) on a 1.7 GHz PC Pentium 4 under Linux. Column 
5 gives the number of terms in the single unsolved equation which in the solution (|70j) below is 
the equation (|7ip. In the remaining 4 columns the number of redundant functions of 6, 5, 4, or 3 
variables is shown. For example, if two functions f(x,y,z) and zg(x) occur always together such 
that a substitution / + zg — > / has the same effect as g — > then g can be set to zero without loss 
of generality. 
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9.3 Conclusions from the test 

The central issue in these runs is, whether integrations (modules 6 and 7) are given a higher priority 
than the formulation of integrability conditions (module 9) or a lower priority. If integrability 
conditions have a higher priority than integrations, as in run 1, then at first a complete differential 
Grobner basis is computed before integrations start. The benefit is that the differential order of 
equations is as low as possible when integrations start (assuming a total degree ordering is used 
in the differential Grobner basis computation). Consequently fewer integrations are necessary and 
fewer functions will be generated which turn out later to be redundant. The disadvantage is that 
the computation of integrability conditions may take very long and blow up the systems size, or 
may even be practically impossible. 

One can attempt to give integrations a higher priority at the price of more redundant functions 
in the solution. This was done in run 2. The benefit may be considerable, only in our small system 
(|69p the Grobner basis computation is not expensive at all, so the advantage of early integrations 
does not become obvious here. But the disadvantage becomes obvious. Integrating higher order 
equations generates more new functions with many of them turning out to be redundant at the end. 

Finally, in the third run we get the best of both previous runs. Here, early integrations use 
syzygies in conservation law form as soon as they become available. The lowered differential order 
of equations reduces the complexity of the remaining Grobner basis computation. Also, because 
with each integration at least 2 equations = P l are satisfied, the number of new functions of 
integration is low and the number of variables these functions depend on is reduced. Consequently, 
only few functions turn out to be redundant in the computed solution as seen in columns 6-9 of 
table 1. 

The following solution is obtained in run 3 after redundant functions have been deleted (by 
hand) leaving 11 functions of 6 variables, 8 functions of 5 variables and 2 functions of 4 variables. 
It is equivalent to the solutions returned in runs 1 and 2. 

c 4 = cioo,x 2 x 3 yiy 2 + ^cioo,^ Xiyt + c i25, x 2 x\y x + ci 25m x 3 yxy 3 + c 133 , X2 x 3 yi + c 133m yiy 3 

+C 2 13,z 3 X l + c 213 m Hi + c 1002/l2/3 - C W9 X 2 X 3 yx + c 170 + c 172 + c 173 x 3 + c 181 y 3 + c 191 (70) 
— C192 — C193X3 — C194 + C200 + C205 — ^ c 65 x 3 y 3 — Cg\X 3 — Cg3 

All functions depend on t, r and in addition on further variables in the following way: 

Cs3(x 2 ,x 3 ,y 1 ,y 2 ), c 8 i(x 2 ,y 1 ,y 2 ,y 3 ), c 173 (x 1 ,x 2 ,y 2 ,y 3 ), c 172 (xi,x 2 ,y 2 ,y 3 ), c 170 (x 1 ,x 2 ,x 3 ,y 2 ), 
ci 94 (xi, x 3 , 2/1,2/3), c 193 (x 1 ,y 1 ,y 2 ,y 3 ), ci 92 (zi, yi, y 2 , y 3 ), c 19 i(x 1 ,x 3 ,yi,y 2 ), c 205 (x 2 ,y 1 ,y 2 ,y 3 ), 
0200(^1,^2, Uu 2/2), cioo(xi,X2,x 3 ), c 125 (x 1 ,x 2 ,y 2 ), c 133 (xi,x 2 ,y 2 ), c 18 i(xi,x 2 ,x 3 ), 
c 2 i 3 (x 2 ,x 3 ,y 3 ), c 23 o(^i,2/i,2/3), c 229 (xi,y u y 3 ), c 2 2 8 (^i, x 3 , 2/1), c 6B (x 2 ,yi), c 109 (x 1 ,y 2 ). 

The function C194 has to satisfy the condition 

= Ci94, X33/1 Xi + 0194,^3 J/1 — 0194,^3 — C228 — C229 — C230^3, (71) 

all other functions are free. The result of the conservation law investigation for the SU(2) chiral 
equation in the form (JBT)|) . (JBTj) is that no other conservation laws with integrating factors of at most 
2 nd order exist. 

More remarks concerning the collaboration of modules: 

• Syzygy based integration can not replace conventional integration. If equations become de- 
coupled then no integrability conditions apply and the equations have to be integrated con- 
ventionally if possible. 
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• The usefulness of conventional integration relies very much on the efficiency of a module for 
the indirect separation (module 8 in the above list). The corresponding implementation in 
Crack will be described elsewhere. 

• The issue of avoiding redundant functions is serious when a system like 1)69)1 is only a subsystem 
of a larger system and the solution of the smaller system is to be substituted in the larger one. 
Redundant functions would complicate the solution of the larger system unnecessarily. On the 
other hand, the identification and deletion of redundant functions using a method described 
in jT3|, is difficult and may be more expensive than the solution/simplification of the system 
itself. This method does not prevent redundancy, it only can identify it in the solution. 

The package Crack is distributed together with the computer algebra system REDUCE. A 
newer version can be down-loaded from http://lie.math.brocku.ca/twolf/crack, 

10 Summary 

An integration method has been proposed that is applicable for linear PDE-systems that admit 
syzygies, i.e. systems which are overdetermined as a whole or contain an overdetermined subsystem. 
It therefore can not replace the straight forward integration of exact PDEs but when applicable it 
has a number of advantages: 

• The information on which the integration is based is taken from syzygies in conservation law 
form. Syzygies are a by-product of the computation of differential Grobner Basis. 

• Because not a single equation is integrated but a number of equations (0 = P l ) at once, fewer 
functions of integration, depending on fewer variables are introduced in the process. 

• The problem of conventional integration to introduce redundant functions when integrating 
with respect to different variables is either prevented or significantly reduced. 

• The new integration produces apart from integrated equations also new syzygies which are 
often the basis for continuing the integration further without having to compute new syzygies 
through a new Grobner basis computation. 

• Syzygy based integration, conventional integration and elimination complement one another 
well in solving overdetermined linear PDE-systems if given the right priorities. 
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Appendix A: Continuation of the introductory example 

In this appendix we continue the introductory example by performing three more syzygy based 
integration steps. The computation is broken up into items. The number(s) at the start of each 
item refer to the line number of the corresponding step in the overview at the end of section El 

(1321) : The remaining system to solve consists of 

= f, xx +f, z (= e 2 ) 

= f,xyz-Cl- (= e 4 ) 
(1331) . (J34|) : satisfying the identity in conservation law form 

= ( e 4;z );CC "I" ( e 2;zj/ e 4)>2 



?5|) : with only 2 derivatives. 
(J36|) : Proceeding as in the first integration step we now identify as the conserved current 

P &&ix f ixxyz QjZ (72) 

P &2ixy 64 f ixxxy Qix (73) 

(I37p : and as the new potential Q we either identify or compute using algorithm DivInt in appendix 
B 

Q fixxy 

+XCi 

(I39|) . (I4()|) : giving the new equation 

= Q - c 2 = f,xxy +xcx - c 2 (=: e s ) (74) 
with the new function of integration c 2 = c 2 (y). 
(I41|) . (I42|) : Equation e 4 is redundant as it turns up purely algebraically in 

= P Q,x = &2ixy 64 C5;x • 

(J43|) : Substitution of e 4 in (|72*j) gives the new identity 

= &2ixxy ~^~^hixx ~^~^hiz • (75) 

(1361) : This is as well a divergence with only two terms 

P X = —G2,xy+G5,x=—f,xyz+Cl = —Q, z (76) 
P z — 65 = f, xxy +XCi — C 2 = Q, x (77) 



(JSU): and the new potential Q 

Q = f,xy +ttCi - xc 2 - 2Cl 



.r 2 



2 
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(J39|) . (J40|) : giving the new equation 

= Q - c 3 = f, xy + —c x - xc 2 - zcx - c 3 (=: e 6 ) (78) 

with the new function of integration c 3 = c^{y). 
(J41|) . (J42|) : Now, equation e$ is redundant as it turns up purely algebraically in 

= P z - Q, x = e 5 - e e , x . 

(J43|) : Substitution of e$ in (|75jl gives the new identity 

= -e2 )a; y +e G , xx +e 6 , 2 . (79) 

(1361) : This is a divergence as well and we will perform the integration cycle one more time with 

P x = -e 2 , y +e 6 , x = -f, yz +xci -c 2 = -Q, z (80) 
P z = e 6 = f\ xy +y Ci ~ XC2 ~ ZCl ~ c 3 = Qix (81) 

(|37j) : and the new potential Q 

Q = f, y +y cx - —c 2 - XZCx + zc 2 - xc 3 
(J39|) . (J40|) : giving the new equation 

= Q - c 4 = /,„ +y c i - y c 2 - ^^Ci + zc 2 - xc 3 - c 4 (=: e 7 ) (82) 

with the new function of integration C4 = c^y). 
(|41|) . (|42|) : Now, equation e6 is redundant as it turns up purely algebraically in P z in (JHTj) 

= P z - Q, x = e 6 - e 7?a; . 

(|43|) : Substitution of e 6 in (|50jl gives the new identity 

= -e 2 , y +e 7 , x . a . +e 7 , 2 . (83) 

The conclusion of this example is shown in section 12.11 below equation (|20j) . As argued there the 
syzygy based integration of equation (J83)) is not advantageous as (J83)) has a conservation law form 
with 3 derivatives instead of two. Instead one rather integrates (|82|) with respect to y and substitutes 
/ in the remaining equation 
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Appendix B: The algorithm DivInt 

The following algorithm computes expressions Q ij (x n JJ) = QM that satisfy D&i = P\ The 
given P % = P l (x™, fj) are assumed to satisfy PjP* = identically in all fj. 



1 Algorithm DivInt 

2 Input variables: x n , functions: f a and conserved current: P l = P i (x n ,fj) 

3 Output Q ij (x n , fj), j >i % satisfying DjQ ij = P\ 

4 E,F % E : list of new additional equations 

5 % F : list of new additional functions 

6 Body % no summation over double indices below 

7 E :— {}, F :— {}, Q ij :— 0, with i,j e 1 p, j >i 

8 

9 % Integrate all terms with functions f a depending on all variables 

10 for i := 1 to (p- 1) do 

11 for j := i + 1 to p do 

12 while P l contains a term a lJ djfJ do % i.e. while any derivative of any f a 

13 % occurs that involves dj 

14 P* -> P i - Dj(a iJ fJ) 

15 P J -> P J ' + D-(a iJ fJ) 

16 q« -> q« + a iJ /y 

17 

18 % Integrate all derivatives involving functions / a not depending on all variables 

19 for % := 2 to p do 

20 for j := 1 to i - 1 do 

21 while P l contains a term a lJ djfJ do % i.e. while any derivative of any f a 

22 % occurs that involves dj 

23 P i -> P* - Dj(a" '/?) 

24 -> - a iJ /y 
25 

26 % Integrate remaining terms 

27 for i : = 1 to p do 

28 if ?V then 

29 % integrate each term a tJ fj of P* with respect to any one x j ^ x % 

30 % preferably one x j with djf a = in the following way: 

31 if d,-/ Q = then 

32 q := fj J a iJ dx j 

33 pi ^ pi- D jq 

34 if j > % then -> + q 

35 else -> - q 

36 else 

37 Introduce a new function /^(x 1 , . . . , x 1-1 , . . . , x p ) 

38 F^FU{f} 

39 P -> E U {0 = 9,-/^ - a iJ fJ} 

40 pi ^ pi - a iJfj 

41 if j > i then Q ij -> Q ij ' + 

42 else Q ji -> - J/ 3 

43 return Q z: >(x n , ff), E (list of new equations), P (list of new functions) 
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Explanation of the algorithm 

Lines 9-16 

This part of the procedure is sufficient if the input expressions P*(x n , fj) do only contain functions 
f a depending on all p independent variables x 1 , . . . , x p . 

A typical example: If an expression P y contains a term f, z then D y P y (no summation) contains 
d y f, z which has to be cancelled by —d z f, y from D Z P Z (no summation) to give = D^P k (summation) 
identically in all fj. This means P z contains —f, y . In this short example the lines 14 - 16 would 
subtract f, z from P y , subtract —f, y from P z and add / to Q yz . There is no principal difference 
between P y containing a term f, z or P y containing a lJ d z fJ. 

As both, P % and P J are updated in lines 14 and 15, j does not run over indices 1 . . . % — 1. Because 
Q n = (Q^ is antisymmetric) there is no need to integrate an i-derivative in P % and therefore j 
starts from i + 1 in line 11. 

If all terms in all P l contain a function f a of all variables then any term in any Q lj occurs 
twice, once with an ^-derivative in P l and once as negative ^-derivative in PK When the program 
completed lines 10 - 16, all P % have the value zero and the solution is found (for i < j, values 
for Qi % follow from the antisymmetry). 

Lines 18 - 42 

The only possibility that after completing lines 10 - 16 not all P l are already zero occurs if some 
f a do not depend on all variables. That is, for example, the case if functions entered the problem 
due to running DivInt previously in earlier integrations. In general, if terms remain in some P l 
which necessarily depend on less than all variables then one can always complete the integrations by 
introducing new functions (collected in a list F in line 38) that have to satisfy additional equations 
(collected in a list E in line 39). In order to minimize the number of additional functions and 
additional equations the lines 19 - 24 integrate terms that are a;- 7 -derivatives in P % (j ^ i) and lines 
31-35 integrate terms by changing the explicit appearance of xK This is shown in the following 
examples. 

Example: Independent variables: x, y, z, initial values: 



containing undetermined functions A, B, C, D, G, H, K, L, M, N, R, S, T, U and W. After complet- 
ing the program up to line 18 the values are 



P 



py 



pz 



A(y, z), y +B(y, z), z +C(y, z) + D(y) + G(z) 
H(x, z), x +K(x, z), z +L(x) + M(x, z) + N(z) 
R(x, y), x +S(x, y), y +T(x) +U(y) + W(x, y) 

Qxy = qxz = Qyz = Q 



P 



py 



pz 



,x 



C(y,z) + D(y)+G(z) 

H(x, z), x +L(x) + M(x, z) + N(z) 

R(x, y), x +S(x, y), y +T(x) + U(y) + W(x, y) 




B(y,z) 
K(x,z). 



After completing the program up to line 26 the values are 



P 



>x 



C(y,z)+D(y) + G(z) 
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py 


= L(x) + M(x,z) + N(z) 


pz 


— t(t) + U(v\ + W(r v) 


Q xy 


= A(y,z)-H(x,z) 


Q xz 


= B(y,z) - R(x,y) 


Q yz 


= K(x,z) - S(x,y). 



The loop beginning in line 27 will integrate the remaining terms in P l . The lines 32 - 35 will 
integrate the terms D, G, L, N, T, U and lines 37 - 42 the terms C, M, W to obtain 

p x = py = p z = o 
Q xy = A(y,z)-H(x,z)+yG(z)-xN(z)+F 1 (y,z) 
Q xz = B(y,z)-R(x,y) + zD(y)-xU(y)-F 3 (x,y) 
Q yz = K(x,z) - S(x,y) + zL(x) -yT(x) + F 2 (x,z) 

with a list F of new additional functions F 1 ^, z), F 2 (x, z), F 3 (x, y) and list E of new additional 
equations 

F\y,z), y = C(y,z) 
F 2 (x,z), z = M(x,z) 
F 3 (x,y), x = W(x,y) 

each in less than 3 variables. 
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